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Abstract. Superconducting quantum circuits are promising systems for 
experiments testing fundamental quantum mechanics on a macroscopic scale 
and for applications in quantum information processing. We report on the 
fabrication and characterization of superconducting flux qubits, readout dc 
SQUIDs, on-chip shunting capacitors, and high-quality coplanar waveguide 
resonators. Furthermore, we discuss the tunability and fundamental symmetry 
aspects inherent to all superconducting qubits, which can be regarded as artificial 
solid-state atoms. Comparing them to their natural counterparts, we discuss 
first and second-order energy shifts due to static control fields. Additionally, we 
present an intuitive derivation of the first- and second-order matrix elements for 
level transitions in the presence of a coherent microwave driving. 
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1. Introduction 

Over the last decade the idea of solid-state based quantum information processing 
has excited scientists and engineers in many disciplines, seeking to combine two 
major developments of the last century, namely quantum mechanics and integrated 
circuit technology [1]. Although the practical challenges to realize solid-state based 
quantum information circuits are huge, the basic building blocks have been successfully 
implemented. Nevertheless, their fabrication and measurement remain demanding 
tasks. At present, most experiments are addressing two key questions. First, how 
can we fabricate, control, and read-out solid-state quantum circuits with sufficiently 
long coherence times, and second, how can we establish a controlled coupling 
between distant circuit parts and transfer the quantum information between them. 
Regarding the first problem, we focus on the fabrication of several key elements of 
superconducting quantum circuits. Their quantum coherent dynamics was discussed 
in our previous work [2, 3, 4]. With respect to the second one, we address symmetry 
aspects of these circuits. In particular, we discuss how symmetry can be broken in a 
controlled way using external control parameters [2]. 
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In this context, superconducting circuits containing nonlinear elements such 
as Josephson junctions or phase-slip centers play an important role. The lateral 
dimensions of these circuits and their constituents can range from 100 nm to a few 
hundreds of jam. Such circuits constitute artificial atoms in the sense that they 
behave similarly to natural atoms in many aspects [5]. Despite the fact that these 
artificial atoms are huge compared to their natural counterparts, they have a discrete 
level structure and exhibit properties unique to the world of quantum mechanics. 
In the simplest case, they form quantum two-level systems, also referred to as 
quantum bits or qubits. Such qubits are particularly attractive for the investigation 
of fundamental quantum phenomena on a macroscopic scale as well as for the use in 
quantum information systems. A major advantage of artificial solid-state atoms over 
natural atoms is their design flexibility and wide tunability by means of external 
parameters such as electric or magnetic fields. Coupling superconducting qubits 
to on-chip microwave resonators has given rise to circuit quantum electrodynamics 
(QED) [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18]. In this prospering field, 
the fundamental interaction between artificial solid-state atoms (qubits) and (single) 
microwave photons is studied. We note that circuit QED is the solid-state analog of 
cavity QED [19, 20, 21] in quantum optics, where natural atoms are coupled to suitable 
three-dimensional cavities. However, in contrast to quantum-optical cavity QED, the 
strong coupling regime, where the qubit-photon interaction coefficient is much larger 
than all relevant decay rates of qubit and resonator, can be achieved easily in circuit 
QED [7, 8]. Moreover, solid-state circuits exhibit an inherent tunability. We have 
recently shown that this can be exploited to break the symmetry of a superconducting 
flux qubit coupled to a lumped-element microwave resonator in a controlled way. 
This allowed us to study multi-photon excitations and to get fundamental insight 
into the interplay of multi-photon processes and symmetries in a qubit-resonator 
system [2]. In another work, we introduced a systematic formalism for two-resonator 
circuit QED, where two on-chip microwave resonators are simultaneously coupled to 
one superconducting qubit [18]. Within this framework, we demonstrated that the 
qubit can function as a quantum switch mediating a tunable interaction between the 
two resonators, which are assumed to be originally independent [18]. 

This article is composed as follows. In section 2, we focus on the fabrication of 
superconducting flux qubits, their readout circuitry and high-quality superconducting 
on-chip resonators. Then, in section 3, we address fundamental symmetry aspects 
of artificial solid-state atoms and compare them to their natural counterparts. As it 
turns out, these results hold for a general quantum two-level system [22]. In particular, 
we spell out the expressions for both flux and charge qubits. 

2. Fabrication of Superconducting Quantum Circuits 

Superconducting quantum circuits are based on flux quantization and Josephson 
tunneling. They have been studied intensively because they have opened a new 
area in fundamental science and because of their potential for the realization of 
quantum information processing systems [5]. Here, we discuss the fabrication process 
of superconducting three-junction flux qubits and the additional circuit elements 
required for readout, coupling and manipulation. 
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2.1. Flux Qubits 

Today, there are three fundamental types of superconducting qubits, namely flux, 
charge and phase qubits [5], depending on the nature of the relevant quantum variable. 
The three-j unction flux qubit [23, 24, 25], we are focussing on here, consists of a 
superconducting loop interrupted by three nm-sized Josephson junctions. Whereas 
two of these junctions with area A and critical current I c are equal, one junction is 
designed to have smaller area a A and critical current al c with a ~ 0.6 — 0.8. In 
an equivalent circuit the Josephson junctions are represented by a normal resistance 
i? n , a junction capacitance Cj and a Josephson inductance Lj = L c /cos(p with 
L c = h/2el c . Here, tp is the phase difference across the junction, h the reduced 
Planck's constant, and e the elementary charge. The characteristic energies are the 
Josephson coupling energy of the junctions, Ej = hl c /2e, associated with the storage 
of a single flux quantum <3>o = h/2e in the Josephson inductor, and the charging 
energy, E c = e 2 /2Cj, associated with the storage of a single electron charge e on the 
junction capacitance Cj. For flux qubits, Ej ^> E c (typically 10 < Ej/E c < 100) 
to ensure that the magnetic flux <E> in the loop is the relevant quantum variable 
and the ground and the first excited state of the qubit have opposite circulating 
currents which produce measurable flux. The minimum temperature for experiments 
on superconducting qubits usually is limited to the base temperature Tb > 10 mK 
of dilution refrigerators. To avoid thermal population of the upper qubit state, the 
level splitting A ge between the ground state |g) and the first excited state |e) should 
be much larger than k^T^, corresponding to an energy and frequency of about 1 fieV 
and 200 MHz, respectively, at Tb = 10 mK. Hence, A ge /h > 2 GHz is desired. Since 
Ag e oc hco p exp(— a^/ Ej / E c ), where co p oc E c ^/Ej/ E c is the plasma frequency and a 
a constant of the order of unity, large E c is required [23, 24]. At the same time one 
has to keep Ej/E c large enough to have well-defined flux states with a measurable 
circulating current. These requirements are demanding regarding junction fabrication 
technology, asking for junctions with high critical current densities J c ~ 10 3 A/cm 2 
and areas down to a few 0.01 jam 2 . Since such small junctions with small parameter 
spread are difficult to fabricate by the well established Nb/AlO^/Nb technology, most 
flux qubits so far are based on A1/A10 X /A1 junctions fabricated by electron beam 
lithography and two-angle shadow evaporation [26]. 

Figure 1(a) shows a scanning electron microscopy (SEM) micrograph of a three- 
junction flux qubit based on Al/AlO^/Al junctions. The qubit is surrounded by the 
readout dc SQUID (Superconducting Quantum Interference Device). Figures 1(d) and 
(e) show an enlarged view of the so-called a-j unction, which has a reduced area a A 
with a = 0.7, and one of the regular qubit junctions, respectively. The area of the 
a-junction is only about 0.02 /mi 2 . The qubits are fabricated on thermally oxidized 
(50 nm SiC^) Si wafers by electron beam lithography using a Philips XL30 SFEG field 
emission SEM and a Raith Elphy Plus nanolithography system. For producing free- 
standing resist masks we used a double-layer resist system. Typical cross-sectional 
views of the resist stencils are shown in Figs. 1(b) and (c). Due to the larger 
sensitivity of the underlay resist, large undercuts can be generated in a controlled 
way. The two-angle shadow evaporation of Al was done by electron beam evaporation 
in a UHV system with a base pressure in the 10 _9 mbar range. The bottom layer 
was thermally oxidized in situ in pure oxygen at p = 2 x 10 -4 mbar. In order 
to achieve high current densities we used small products L of oxygen pressure and 
oxidation time ranging between 0.25 and 0.26 mbar • s. With these values we achieved 




200 nm 



Figure 1. (a) SEM micrograph of a three-junction flux qubit (inner loop) 
surrounded by the readout dc SQUID (outer loop with two junctions). The 
white rectangles mark the regions of the a-} unction (upper) with a reduced 
area (a = 0.7) and a regular junction (lower) shown on an enlarged scale in 
(d) and (e), respectively, (b) and (c) show cross-sectional views of double- 
layer resist structures obtained with different exposure times. Large undercuts 
can be obtained, resulting in free-standing resist parts required for the shadow 
evaporation technique. 



J c (50mK) c± 1500 ± 500 A/cm 2 with a high yield of 85%. By optimizing the electron 
beam lithography process we could reduce the spread in the junction area down to 
10%. We note that we can increase or decrease J c by decreasing or increasing L, 
respectively. For example, by increasing L to above 10 mbar • s the J c values can be 
reduced to below 10 A/cm 2 . 

Figure 2 shows the current- voltage characteristics (IVC) and the flux dependence 
of the critical current of a readout dc SQUID fabricated with the process described 
above. From I c = 1.44 /iA and the measured junction area A = 0.035 jam 2 we 
obtain J c = 2kA/cm 2 . Furthermore, from R n = 169 Q the product I c R n = 243 /iV 
is obtained. The measured gap voltage V g ~ 360 /aV is close to the BCS value 
V g = 2A/e = 3.53/cBT c /e = 365 /iV using T c = 1.2 K for the critical temperature 
of Al. Furthermore, the measured I c R n product agrees well with the Ambegaokar- 
Baratoff value [27] 7rV^/4 = 287/iV. Using the specific junction capacitance C s = 
100 ±25 fF /jam 2 [28], the capacitance of the qubit junctions with area A = 0.02 jam 2 is 
estimated to be Cj ~ 2fF corresponding to E c = 40 /aeV or, equivalently, T = 0.46 K. 
The subgap structures in the IVC in Figure 2(a) most likely originate from LC 
resonances as has been discussed in more detail recently [28]. 

We note that the fabrication of sub-micron sized Josephson junctions with two- 
angle shadow evaporation is very common since it is easy to use. However, there 
are also drawbacks. In particular, the reproducibility and parameter spread of the 
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Figure 2. (a) Current-voltage characteristics of a readout dc SQUID based on 
A1/A10 X /A1 tunnel junctions with A = 0.035 /im 2 at T = 50 mK. (b) Critical 
current I c of a readout dc SQUID as a function of the applied magnetic flux <E>* X 
threading the SQUID loop at T = 540 mK. The inset shows an enlarged view 
of the region around <E» e x = — 7<E>o/2, where the qubit persistent current changes 
sign. Here, 7 = 1.37 is the ratio of the loop areas of the SQUID and the qubit. 



junctions will be finally limited by the non-planar junction geometry, resulting in 
a not well defined junction area. The related fluctuations in the critical current 
become an increasing problem when one is going to more complex quantum circuits. 
However, more severe is critical current noise due to impurities and trapping sites 
in the tunneling barrier leading to decoherence. Therefore, in the long run a planar 
junction geometry and the use of epitaxial electrode and barrier layers may be the 
optimum choice. 

The J C (<I>* X ) dependence of a readout dc SQUID is shown in Fig. 2(b). The 
measured curve is close to the ideal | cos(7r<I>£ x /<I>o)| dependence expected for a dc 
SQUID with /?l = 2ttL s I c /^q = L s /L c <C 1, for which the Josephson inductance 
L c dominates the geometric inductance L s . However, there are small deviations, 
originating from the additional flux due to the persistent current I p circulating clock- 
or counterclockwise in the qubit loop. In the region close to the applied flux <l> ex 
corresponding to $o/2 in the qubit loop, 7 p changes sign. This results in a step-like 
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Figure 3. (a) Normalized switching current J sw of the readout dc SQUID plotted 
versus the frustration / = ^ex/^o in the region around the qubit step at T = 
50 mK. The qubit is irradiated by a microwave signal (y = 8 GHz). The dashed 
curve shows the theoretically expected curve without microwave irradiation for 
Ag e /h = 5.2 GHz. The inset shows the switching current distribution in a color 
coded image, (b) Qubit transition frequency E ge /h (symbols) plotted versus the 
frustration /. The line is a numerical fit to the data yielding A ge /h = 5.2 GHz 
and J p = 450 nA. 



feature superimposed on the regular I c ($* x ) dependence. This so-called qubit step is 
shown in the inset of Fig. 2(b). With the ratio of 1.37 between the SQUID and the 
qubit loop we estimate <£ ex = — 0.685 <3>o m good agreement with the experiment. 

The qubit level structure has been investigated by microwave spectroscopy. 
Figure 3(a) shows the normalized switching current of the readout SQUID around 
the qubit step plotted versus the frustration / = ^ ex /^>o- Here, <£> ex is the applied 
flux threading the qubit loop. Note that the switching of the SQUID from the zero- 
voltage into the voltage state is by quantum tunneling and hence is a statistical 
process. Therefore, the normalized switching current I S w(f) is obtained from the 
switching current distribution measured for every /. This distribution is shown in a 
color-coded image in the inset of Fig. 3(a). The J S w(/) curve shows a pronounced 
peak-dip structure due to resonant transitions between the qubit levels induced by 
the applied microwave irradiation with frequency v = 8 GHz. Note that multi-photon 
transitions are observed depending on the applied microwave power. The position of 
the peaks/dips (one-photon transitions) as a function of frequency and frustration is 
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Figure 4. (a) Optical micrograph of the readout and bias circuit of the dc SQUID 
and the end of the coplanar waveguide which is shorted to ground to form an 
antenna. In (b) an enlarged view of the region around the qubit and SQUID loop 
with the shunting capacitors is shown. The equivalent circuit is depicted in (c). 



shown in Fig. 3(b). The data follow the expected E ge = ye 2 + A 2 e dependence, where 

e = 2ip<3>o (7+0.5) is the energy bias and A ge the level splitting at the degeneracy point 
/ = —0.5. Fitting the data yields J p = 450 nA and a qubit gap of A ge /h = 5.2 GHz. 

2.2. Bias and Readout Circuit 

For the manipulation and readout of flux qubits as well as for the attenuation of 
external electrical and magnetic noise additional circuit elements such as resistors, 
capacitors, inductors, microwave lines, antennas etc. are required. Some of them have 
to be fabricated on-chip. For the fabrication of these usually larger circuit parts we 
used a mix&match process, where high resolution electron beam lithography patterns 
are precisely aligned into existing patterns made by optical lithography. For the 
alignment markers, the leads and the coplanar microwave lines we used 20 nm thick Au 
layers. These layers are fabricated by optical lithography, electron beam evaporation 
and a final lift-off process prior to the e-beam process for the qubits described above. 
During the subsequent Al evaporation process, parts of the Au pattern are covered by 
Al. The small thickness of the Au layer allows to avoid the interruption of the Al film 
at the edges of the Au layer. 

Figures 4(a) and (b) show optical micrographs of the circuit elements used for 
biasing the readout dc SQUID and measuring the SQUID voltage. As shown by the 
equivalent circuit in Fig. 4(c), the SQUID is shunted to ground by on-chip Al/AlO^/Al 
capacitors with the total capacitance C. In combination with the on-chip resistors 
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in the dc SQUID bias and voltage lines, which are realized by long Au lines, they 
form a low-pass filter constituting the main component of the qubit electromagnetic 
environment. Details on the shaping of the electromagnetic environment as well as 
the effect of a resistive and capacitive coupling of the dc SQUID are given elsewhere 
[4, 3]. The ground plane of the Al/AlO^/Al capacitors is deposited by electron beam 
evaporation prior to the qubit fabrication process. It is thermally oxidized in ambient 
atmosphere for 3 hours at 100° C to achieve a thick oxide layer with sufficiently small 
leakage current. Specific capacitance values between 12 and 15fF//im 2 and resistance 
times area products larger than 100 GQjam 2 could be achieved. 

It is evident from Fig. 4(c) that the two capacitors with total capacitance C/2 
together with the two inductors with total inductance 2L form an LC resonator 
coupled inductively to the qubit loop. The typical resonance frequency v x = uo T /2tt 
of the LC resonator is in the GHz range and is comparable to the qubit transition 
frequency. Furthermore, the coupling strength g/2ir is of the order of 100 MHz [2, 11]. 
Therefore, the coupled qubit-resonator system allows to study fundamental questions 
of the flourishing field of circuit QED [7, 8, 11, 12, 13, 14, 15, 29]. Here, the interaction 
of solid-state quantum systems with single microwave photons is of particular interest. 
Recently, we have used a circuit very similar to that shown in Fig. 4 to observe one 
of the key signatures of a two-photon driven Jaynes-Cummings model, which unveils 
the upconversion dynamics of a superconducting flux qubit coupled to an on-chip 
lumped-element microwave resonator [2]. 

2.3. Microwave Resonators 

The quality factor Ql of the LC resonators discussed above is only of the order of 
100. Therefore, the decay rate n = is t /Ql of the microwave photons from the resonator 
is large. In order to access the strong coupling regime in circuit QED experiments, 
g/27r ^> k is required. With gj2i\ ~ 100 MHz and uj y /2tt ~ 10 GHz, microwave 
resonators with Ql ^ 10 4 have to be designed. Such resonators can be realized by 
simple transmission line geometries [30]. 

We have fabricated coplanar waveguide (CPW) resonators based on 200 nm thick 
Nb films on thermally oxidized Si and sapphire substrates. An optical micrograph 
of such a resonator is shown in Fig. 5(a). The Nb films have been deposited by dc 
magnetron sputtering in Ar. The patterning was done by optical lithography and 
reactive ion etching using an Ar/SF6 mixture. The center conductor of the CPW 
resonators is 20 jam wide, separated from the lateral ground planes by a 12 jam wide 
gap to obtain a wave impedance of 50^. The resonator is coupled at both ends to 
a coplanar transmission line by interdigital or gap capacitors shown in the insets of 
Fig. 5(a). The strength of the capacitive coupling determines the loaded quality factor 
Q^ 1 = Q 1 + QeAi where Q is the internal quality factor of the uncoupled resonator 
and Qext is determined by the value of the coupling capacitors C ext and the load 
impedance (50 Q in our case). By proper choice of C ex t the loaded quality factor can 
be designed to an optimum value satisfying on the one hand the condition g ^> k and 
allowing on the other hand an as high as possible measuring rate (ex n). The loaded 
quality factor Ql has been extracted from S21 measurements as shown in Fig. 5(d). In 
Fig. 5(c) we have plotted Ql as a function of the coupling capacitance C ex t, calculated 
using FastCap®. The measured dependence nicely follows the theoretical prediction 
[31] showing that Ql can be deliberately designed over a wide range. As shown in 
Fig. 5(d), internal quality factors Q > 10 5 can be achieved with the niobium CPW 
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Figure 5. (a) optical micrograph of a CPW microwave resonator. The 
insets show different coupling capacitors located at the positions marked by 
the rectangles, (b) Picture of the gold-plated copper sample box containing a 
resonator chip. The rectangles mark the positions of the coupling capacitors, 
(c) Measured Ql (symbols) plotted versus the coupling capacitance for a 6 GHz 
niobium CPW resonator at T = 1.5 K. The dashed line shows the theoretical 
prediction [31]. (d) Transmission spectrum of an undercoupled 1.25 GHz niobium 
CPW resonator measured at 1.5 and 4.2 K. 



resonators already at 1.5 K. 



3. Tunability, Coupling Strength and Symmetry Breaking in 
Superconducting Circuit QED 

The field of superconducting circuit QED has attracted considerable attention over 
the last years, since it allows to study the interaction of microwave photons and solid- 
state quantum circuits on a quantum level. In the following we consider a classically 
driven quantum system [32] consisting of a qubit coupled to a harmonic oscillator. 
We note that the subsequent discussion is quite general. It is not restricted to the 
flux qubits discussed in section 2, but is valid for a general quantum two- level system 
including both superconducting flux and charge qubits. The first-order hamiltonian 
of this system in the diabatic basis of the qubit eigenstates {| |), | |)} is given by [2, 9] 

H = e -d z - ^y°x + {a} a + 

+ hg a z (a^ + a) H a z cos cot + rj (a^ + a) cos out . (1) 
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Here, the Hamiltonian H q = (e/2)a z — (A ge /2)a x describes the qubit and H T = 
fyjO T (a)a + 1/2) represents a quantum harmonic oscillator with the photon number 
states |0), |1), 1 2), . . .. a) and a are the boson creation and annihilation operators and 
a x and a z the Pauli operators. The level splitting at the degeneracy point A ge > by 
definition. The interaction term between qubit and resonator is H^ T = hgd z (a) + a), 
and i^ m ,q = {£l/2)a z coscut and H mj = rj(aJ + a)cosut describe the qubit and the 
resonator microwave driving, respectively. Furthermore, uj v /2it is the fundamental 
resonator frequency, g/2ii the qubit-resonator coupling strength, and ft and n are the 
qubit and the resonator driving strength. At the qubit optimal point and, under some 
approximations, also away from it, the hamiltonian (1) is equivalent to that describing 
a driven system consisting of an atom coupled to the light field in an optical cavity. 
Therefore, solid-state circuit QED can be considered the circuit analogue of cavity 
QED. We note that the hamiltonian (1) can be applied to both flux and charge qubits 
by choosing an e = 2/ p £<I> ex and e = 2V p 5Q ex , respectively. Here, V v = e/C g is 
the voltage associated with an elementary charge on the gate capacitance C g , 5& ex 
and SQ ex are the applied gate flux and charge relative to the degeneracy point. The 
latter is located at half a flux quantum and half a Cooper pair charge, respectively. In 
general, one can switch between the expressions for flux and charge qubits by replacing 
the quantities current, flux and inductance by voltage, charge and capacitance. 

Most interesting in circuit QED is the so-called strong coupling regime, g/2ir ^> 
ft, 7, where n is the resonator decay rate and 7 is the largest of the qubit decay rates. In 
this situation, the reversible exchange of microwave photons between a solid-state two- 
level quantum system (qubit) and a resonator represents a coherent oscillatory process. 
Regarding quantum information processing, this means that quantum information can 
be exchanged back and forth many times. To achieve the strong coupling regime, 
the qubit-resonator coupling strength g/2ir must be of the order of 100 MHz because 
presently 7 > 1MHz, i.e., the decay rate of the qubit into other channels is still 
large for superconducting qubits. As demonstrated in section 2, the fabrication of 
suitable resonators with n < 1 MHz is possible. Compared to optical cavity QED, 
such large values for g/2ir can be obtained more easily in superconducting circuit 
QED for two reasons. First, the electric and magnetic dipole moments, ji e \ and 
Mmag, of superconducting charge and flux qubits, respectively, are much larger than 
those of natural atoms. Second, the mode volume V m of transmission line microwave 
resonators (cf. Fig. 5) can be made much smaller than for 3D cavities in quantum 
optics, resulting in large electric or magnetic fields. Note that the fields associated 
with vacuum fluctuations are Eq = V / huj r / / 2eoKn and Bq = y/ /HQhcj r /2V m . Since 
mode volumes as small as 10 _12 m 3 seem to be feasible, the electric and magnetic 
fields associated with vacuum fluctuations are of the order of 0.1 V/m and 1 nT for 
uj v /2it of a few GHz. 

In the following we first address the tunability of the energy levels of 
superconducting charge and flux qubits by static electric and magnetic fields. The 
derived energy shifts are equivalent to the dc Stark and Zeeman shifts in natural 
atoms. Then, we discuss the atom-resonator coupling strength and some fundamental 
symmetry aspects of superconducting circuit QED. In our discussion we use a 
semiclassical approach, treating the microwave driving as classical fields. 
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3.1. Tunability of Charge and Flux Qubits 

In this section we compare artificial superconducting atoms (charge and flux qubits) 
to natural atoms (we use hydrogen for simplicity) regarding their tunability by static 
fields. In Table 1 we have summarized the first- and second-order energy shifts 
AEgJ^g and AE^ m due to electric and magnetic fields, SE ex and 5B ex , applied 
via suitable control gates. Discussing the first-order terms, we consider the charge 
qubit at the degeneracy point as a system of two degenerate levels | j) and | |) 
corresponding to adjacent Cooper-pair number states with Cooper pair numbers N 
and N + 1. At this point, the polarization charge Q ex = 2e7V ex induced by an 
applied gate voltage corresponds to 2e(N + 1/2) for both charge basis states. We 
define SN ex — 7V ex — (TV + 1/2) as the deviation from the Cooper pair number at 
the degeneracy point. Therefore, at the degeneracy point we have 5N ex = and, 
correspondingly, 5E ex = 0, where SE ex = 5Q ex /C g d = 2e8N ex /C g d is the electric field 
applied via a control gate, which just shifts the Cooper pair number on the island 
by 5N ex away from the degeneracy point. Here, C g is the gate capacitor and d the 
distance between the control gate and the island. Our discussion shows that we can 
consider the charge qubit at the degeneracy point as a degenerate two-level system 
in a zero effective control electric field SE ex = 0. The finite coupling of the charge 
states | |) and | j) results in the symmetric and antisymmetric superposition states 
|g) = (ll) + lT))/v / 2 and |e) = (| |) - | T))A/2> respectively, with an energy splitting 
Ag e of their eigenenergies. The situation for the flux qubit is completely analogous. 
At the degeneracy point it can be considered as a system of two degenerate levels | j) 
and | j) corresponding to the persistent currents ±J p circulating in the qubit loop. 
The polarization flux <l> ex = 7V ex ^> is generated by a magnetic field B ex applied via a 
control gate. Again, at the degeneracy point we have SN ex — and SB ex = 0. Here, 
5B ex = 5$> ex /A = 6N e ^o/A is the magnetic control field, which just changes the flux 
number in the loop of area A by SN ex . Evidently, at the degeneracy point we can 
consider the flux qubit as a two-level system in a zero effective magnetic control field 
5B ex = 0. As for the charge qubit, the finite coupling of the flux states results in 
superposition states with eigenenergies separated by A ge . 

We now can discuss the interaction energy of the charge and flux qubit with 
static control fields SF ex — {5E ex ,SB ex }. Since both the electric and magnetic dipole 
moments in the two basis states | j) and | j) have opposite sign, the interaction 
operator in the basis {| j), I T)} can be written as 

W = -^5F ex a z (2) 

with the dipole moments /i = {/i e i,/i mag } given by 

Mei = e d (3) 

Mmag = IpA . (4) 

We immediately see that the first-order energy shifts of the uncoupled levels | j) and 
| |) are given by the products 

A^ 1} = ± /i el • SEe* (5) 

A^ag = ± AWg • 5Bex • (6) 

For natural atoms, the corresponding energy shifts are the linear Stark and Zeeman 
shifts. Note that natural atoms and molecules with an inversion center in a non- 
degenerate electronic state do not have a permanent electric dipole moment and 
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Table 1. Characteristic first- and second-order shifts of the energy levels 
due to external electric and magnetic fields SE ex and (5£> ex for natural atoms 
(hydrogen) compared to those of superconducting charge and flux qubits. For 
the estimates of the electric and magnetic dipole moments fj, e \ and fi ma g as well 
as the polarizabilities Xe\ an d Xmag we used the typical values d = 0.1 — 1 jum 
for the distance between island and gate electrode (charge qubit), as well as 
A = 5 — 10 /im 2 and I p = 100 — 500 nA for the loop area and the persistent 
current (flux qubit), respectively. Furthermore, /io = 1-26 x 10 _6 Vs/Am is 
the vacuum permeability, eo — 8.85 x 10 — 12 As/Vm the vacuum permittivity, 
a B = 5.29 x 10 _11 m the Bohr radius, m = 9.109 x 10 _31 kg the electron mass, 
and e = 1.602 x 10 _19 As the electron charge. 

energy natural atom charge qubit flux qubit 

shift 

AE^ ±/i e l • SE ex ±/i e l • ^ex 

(Mel) =0 Mel = ed 

~ 10- 26 -10- 25 Cm 

A£/mag i/^mag ' ^ex iMmag ' ^ex 

Mmag — MB Mmag = -^p^- 

- 10 4 -10 5 mb 

Xei - 47re a| Xei - ~ 2 ^^ 

~ 10- 41 Cm 2 /V ~ 10- 28 -10- 26 Cm 2 /V 

a 77>( 2 ) AR 2 AR 2 

a -^mag 2^ ma S ' aD ex 2^ ma S ' aD ex 

^ e 2 a| ^ (/ p A) 2 

Xmag — 6m Xmag — ^ A ge 

- IO" 29 Am 2 /T - 10- 13 -10- n Am 2 /T 



therefore have ^ = 0, that is, they do not show a linear Stark effect. As indicated 
by the numbers listed in Table 1, the dipole moments of the qubits are by many orders 
of magnitudes larger than those of natural atoms. For the charge qubit, fi e \ = ed is 
determined by the distance d between the island and the gate electrode and can be as 
large as 10 4 to 10 5 debye (ldebye = 3.3 x 10 _30 Cm). In other words, it exceeds the 
typical values for natural atoms by 4 to 5 orders of magnitude. The same is true for 
the magnetic dipole moment, which is given by the persistent current I p circulating 
around the qubit loop and the loop area A. It typically amounts to 10 4 to 10 5 Bohr's 
magnetons compared to only a single Bohr's magneton for natural atoms. Obviously, 
the large dipole moments of the qubits are caused by their much larger physical size. 
Note that d and A are of the order of /im and /im 2 for the qubits, but only A and 
A 2 for natural atoms. Of course, the large dipole moments of qubits also have a 
disadvantage. They make the qubits more sensitive to electric and magnetic noise. 
We also note that large dipole moments are obtained for Rydberg atoms [33, 34, 35], 
which are highly excited atomic states with much larger atomic diameter. 
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So far we have considered only the uncoupled basis states. The finite coupling of 
| |) and | |) results in the superposition states [22] 

\ g )= +C os^e + ^ 2 U)-sin^e-^ 2 |T) (7) 

\e) = + sin 6 -e +i ^ 2 \ |) + cos -e~^ 2 \ |) (8) 
with the level spacing 

+ A 2 e . (9) 
Here, tan# = A ge /e and e = 2/i £F ex is the energy bias given by 

(2e) 2 

e e i = 2ed ££ ex = 2J7V ex (10) 

e mag = 2/ p A ££ ex = H 2J7V ex (11) 

for the charge and flux qubit, respectively. The large dipole moments of the qubits 
allow to considerably tune E ge by applied control fields. Already fields of the order 
of 10 3 V/m and 10 _4 T are sufficient to cause changes of the order of A ge , whereas 
for natural atoms the equivalent shifts are usually negligibly small. It is important 
to note that the expectation values are (g\a z \g) = —cosO and (e\a z \e) = cosO. Since 
cos# = at the degeneracy point, the expectation value of the dipole moments vanish 
for the superposition states. Therefore, the energy shifts AE^ and AE'mag disappear, 
making the charge and flux qubit insensitive to low-frequency fluctuations (e.g. 1/f 
noise) of the electric and magnetic field at this operation point, respectively [3, 4, 36]. 
Nevertheless, there are significant second-order contributions left, originating from the 
large polarizabilities of the qubits as discussed in the following. 

The second-order terms can be expressed by using the polarizabilities of the 
natural and artificial atoms. For natural atoms with spherical symmetry, the 

(2) 1 

polarizability tensor is isotropic, giving the quadratic terms AE^ } J = — ^Xei * 
5E* X and AE^l g = 

— ^Xmag ' ^ex5 where Xei an d Xmag are the electric and 
magnetic polarizabilities, respectively. Using second-order perturbation theory, the 
polarizabilities of the qubits can be written as 

(^|fel,mag|e)(e|/iel,mag^) 9 x 
Ael,mag z rp V iZ / 

^ge 

with E ge = E e — E g1 /i e i = edd zi and /2 mag = I p Aa z . Because (g\a z \e) = (e\a z \g) = 
sin#, we obtain the polarizabilities 



and the second-order energy shifts 
,(2) _ (edf 



.2 , 



A^r =^ < as) 



(2) {I p Af rin'g 

mag 771 ex • V xv v 



£ge 
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Interestingly, the second-order contributions are maximum at the degeneracy point, 
where the first-order contributions vanish, and decrease moving away from it because 
of the sin 2 term. Close to the degeneracy point we can use the approximations 
sin# ~ 1 and E ge c± A ge and obtain the polarizabilities Xei — — 2(e<f) 2 /A ge and 
Xmag — — 2 (ip A) 2 /Ag e , which are listed in Table 1. Again, due to the much larger size 
of the qubits the polarizabilities are by many orders of magnitudes larger than those of 
natural atoms. Since in electrostatics the electric and magnetic energy contributions 
also can be expressed as ^CqcPSE^ and \{A 2 /Lq)5B' 2 x , the polarizabilities can be 
associated with a quantum capacitance Cq — 2e 2 /A ge and a quantum inductance 
L Q = A ge /2/ 2 , respectively [37, 38, 39, 40]. 

3.2. Coupling Strength and Symmetry Breaking in Circuit QED 

In the previous subsection we have seen that charge and flux qubits show a large 
tunability of their energy levels due to their large dipole moments and polarizabilities, 
allowing to tune the level separation over a wide range by static electric and magnetic 
gate fields. In the following we address the coupling strength with dynamic electric 
and magnetic fields as well as some fundamental symmetry aspects of superconducting 
circuit QED. In our analysis we do not only focus on flux qubits discussed in section 2.1. 
We rather show that the arguments given in the following are quite general and can 
be straightforwardly applied also to charge qubits. 

It is well known [24] that the potential of the three-junction flux qubit can be 
reduced to the one-dimensional double- well U(&) shown in Fig. 6(a). Here, $ is the 
total magnetic flux in the qubit loop. It originates from the externally applied flux <3> e x 
and the flux due to the current circulating in the loop, having an expectation value 
I p (a z ). We first discuss the situation at the degeneracy point, where the external flux 
$ ex equals = <&o(n + 1/2). Here, n is an integer and we only consider n = 
in the following. Then, the two wells are symmetrical with respect to <!> = $o/2 and 
I p ((Tz) vanishes. In Fig. 6(a) we also plot the eigenfunctions |g) and |e), which, at 
the degeneracy point, are fully symmetric and antisymmetric superpositions of the 
two basis states | j) and | j) with clock- and counterclockwise circulating persistent 
currents ±/ p . For a wider region around the degeneracy point, the energy difference 
£ge(e) (cf. eq.(9)) between the two levels is shown in Fig. 6(b). At the degeneracy 
point, the energy bias e = 2I p ASB ex = 2I p 6^ ex is zero resulting in ^g e = A ge . The 
probability to observe | j) (| |)) is just 0.5 at the degeneracy point and approaches 1 
or (0 or 1) when decreasing or increasing the applied flux. For the charge qubit the 
situation is completely analogous. At the degeneracy point, the ground state |g) and 
the excited state |e) are fully symmetric and antisymmetric superpositions of the two 
basis states | j) and | j) corresponding to neighboring Cooper pair number states on 
the superconducting island. 

Considering symmetries, we immediately see from Fig. 6(a) that the potential 
U of the flux qubit is fully symmetric with respect to the variable 5& = $ — &o/2 
for an external flux bias <£ e x = 3>o/2. That is, at the degeneracy point we have 
U(5$>) = U(— S$>). We note that 5$ is associated with the phase differences 71 and 
72 across the two larger junctions as (72 — 7i)/2 = 2tt(S^/^o)- In complete analogy, 
U(SQ) = U(—SQ) for the charge qubit at an external charge bias of Q ex = 2e/2. 
Here, 5Q = Q — 2e/2 is the quantum variable conjugate to 5<f>. In both cases the qubit 
eigenstates |g) and |e) have even and odd parity because they are fully symmetric and 
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Figure 6. (a) Sketch of the double- well potential U of a flux qubit plotted versus 
the total flux <E> in the qubit loop for an externally applied flux <E> e x = ^o/2 
(solid green line). Also shown are the two eigenfunctions of the ground state 
(symmetric, red) and the first excited state (antisymmetric, blue). The dotted 
and dashed lines show the qubit potential for an external flux larger and smaller 
than <E>o/2. (b) The energy E of the two superposition states |g) = l/y/2(\ J,) + | f)) 
and |e) = l/y/2(\ I) — \ T)) versus the energy bias e. The dashed lines show the 
energies of the uncoupled basis states. 



antisymmetric superpositions of the basis states. Then we can write 

n |g> = + |g> (it) 

fi |e> = - |e> , (18) 

where II is the parity operator. Moving away from the degeneracy point, the potential 
is tilted to the left or right as indicated by the dotted and dashed lines in Fig. 6(a). 
Then, symmetry is no longer well defined and |g) and |e) are no longer fully symmetric 
and antisymmetric superpositions of the basis states. In this way, we we can break the 
symmetry of the flux qubit in a controlled way using an external control parameter, 
namely the static magnetic flux 5$> ex = <I> ex — <I>o/2 generated by a magnetic field 
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applied via a control gate. For charge qubits, correspondingly the control parameter 
is the excess charge 5Q ex = 2e5N ex = 2e(iV ex — 1/2) induced on the island by an 
electric field applied via a gate capacitor. 

The controlled symmetry breaking by static applied gate fields has far-reaching 
consequences on the selection rules for level transitions. At the degeneracy point, 
these selection rules can be intuitively understood by the fact that the qubit potential 
exhibits mirror symmetry, U(5<$>) = U(—5&) and U(SQ) = U(—5Q) for the flux 
and charge qubit, respectively [cf. Fig. 6(a)]. Hence, the interaction operator of 
the one-photon driving initiating level transitions between the fully symmetric and 
antisymmetric states \g) and \e) must be odd with respect to the variable S$ ex 
or 5Q ex [2, 41]. Such odd-parity interaction operators are the dipole operators 
[cf. eqs.(2), (3), and (4)] for the flux and charge qubit, respectively, because 

^2 g ,e K)(^lMmag,ei|j)0'| cx d x and, hence, the anticommutator {II, /2 magje i} = 0- In 
contrast the interaction operator for the two-photon driving is even [2]. Hence, if the 
interaction operator is a dipole operator, at the degeneracy point only one-photon 
transitions are allowed, whereas two-photon transitions are strictly forbidden. These 
selection rules correspond to those for electric dipole transitions in natural atoms 
[22]. Away from the degeneracy point, the qubit potential no longer exhibits mirror 
symmetry with respect to 5<f> and SQ for the flux and charge qubit, respectively, and 
the superposition states \g) and \e) no longer are purely symmetric and antisymmetric 
states. Then, the strict selection rules valid at the degeneracy point do no longer 
apply and one- and two-photon transitions can coexist. 

To gain a better insight into this issue one has to explicitly discuss the matrix 
elements for level transitions in a flux qubit induced by an external microwave 
driving SB ex (t) = S Bq cos cot and SE ex (t) = SEq cos cot of the flux and charge qubit, 
respectively. In this context, we first note that for the flux qubit the transitions 
induced by 5B ex (t) are of electric dipole type. This can be understood by recalling 
that for the flux qubit the magnetic flux is equivalent to the spatial coordinate, 
whereas Q is equivalent to the momentum. In contrast, for the charge qubit, the 
charge Q is equivalent to the spatial coordinate and the flux $ to the momentum. 
With the dipole operator D(t) = ^ cos cot a z and using (g\a z \e) = sin# we obtain the 
matrix element for one-photon transitions [22] 

(g\D(t)\e) = ^ (g\S z \e) coscot 

= — sin# cos cot , (19) 

where we define Q = 2I p ASBq and ft = 2V p CdSEo = 2ed5Eo for the flux and 
charge qubit, respectively. Hence, the modulus of the matrix element for one-photon 
transitions is given by 

D« = |sin0. (20) 

For both the flux and charge qubit the transition matrix element is proportional to 
sin# = A ge /E ge . Hence, it is maximum at the degeneracy point (sin# = 1) and 
continuously decreases when moving away from it. 

The treatment of the two-photon transitions is more complicated. For 2co 
co ge = E ge /h (resonant two-photon driving) the level transitions can be described as 
a two-step process and a detuning of E ge /2 [42]. Then, the matrix element for the 
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Figure 7. Selection rules for one- and two-photon transitions in a flux qubit 
having a potential with (left) and without (right) mirror symmetry. In the 
symmetric situation (left) the two-photon process is forbidden. Since two-photon 
transitions are also forbidden for the harmonic oscillator, qubit and resonator can 
be excited only resonantly at frequency uj\ = uj ge = uj v . In contrast, in the case of 
broken symmetry (right) the qubit can be excited both by one- and two-photon 
processes. Therefore, also for u>2 = <^ ge /2 (two-photon process) the excited qubit 
state \e) can be populated. Due to the qubit-resonator coupling g/2it a state 
transfer to the resonator is achieved. The resonator state, in turn, decays by 
emission of a photon of twice the frequency (frequency up-conversion). 



two-photon transition is obtained to 

(g\D(t)\e)(e\D(t)\e) = Q 2 (g\a z \e)(e\a z \e) 2 

E ge /2 4 E ge /2 ' [ } 

With (g\a z \e) — sin#, (e|<r z |e) = cos#, cos 2 ut = ^(1 + cos2o;t) and E ge = A ge / sin# 
we obtain 

(g\D(t)\e)(e\D(t)\e) . 2 n n . 

x — v J ,1 = -r-r— sin 2 cos (9(1 + cos 2ut) . (22) 

E ge /2 4 A ge 

From this expression, we identity the modulus of the matrix element for two-photon 
transitions 

£>(2) = sin 2 cos 0. (23) 

4A ge 

The transition matrix elements for both flux and charge qubits are now proportional 
to sin 2 cos = A ge e/E ge . Because e = at the degeneracy point, no two-photon 
transitions are allowed there in agreement with the qualitative discussion above. 
However, moving away from the degeneracy point the matrix element becomes 
finite, making a two-photon driving of the qubit possible. This has been demonstrated 
in a recent experiment [2], where a flux qubit coupled to an LC resonator has been 
studied. 
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The possibility to deliberately tune the level spacing and the symmetry properties 
of flux and charge qubits by an external control field can be used to study the effect 
of symmetry breaking in superconducting circuit QED systems. The key effect of 
symmetry breaking on selection rules and matrix elements for one- and two-photon 
transitions in a coupled flux qubit-resonator system is sketched in Fig. 7 for a three- 
junction flux qubit. The resonator, a harmonic oscillator, can never accept a two- 
photon driving. Furthermore, for an external flux bias <£ ex = <l>o/2 the qubit potential 
has perfect mirror symmetry (left) and the selection rules discussed above hold. Then, 
qubit and resonator can be excited only resonantly at frequency uj\ = co ge . By 
changing the flux bias away form <3>o/2, the mirror symmetry of the potential is broken 
(right). In this case, also for two-photon driving (c^2 = ^ge/2) the excited qubit state 
\e) can be populated. The strong qubit-resonator coupling g/2ii results in a state 
transfer to the resonator, which subsequently decays by emission of a photon of twice 
the frequency. This frequency up-conversion mechanism for microwave photons in 
a coupled qubit-resonator system has been studied in a recent experiment [2]. We 
note that the occurrence of this process is a direct consequence of the underlying 
symmetry properties of the system, which can be chosen on purpose by an external 
control parameter (magnetic flux in our case). The symmetry properties of the qubit- 
resonator system have an interesting consequence. Depending on the choice of <3> ex 
and the microwave frequency o;/27r, one can choose to drive either only the qubit, 
only the resonator, or both the qubit and the resonator. This selective driving once 
more demonstrates the benefits of the tunability of solid-state circuit QED systems. 
Besides its role in fundamental research, the controlled symmetry breaking in artificial 
quantum systems has good prospects of use in numerous applications. Particular 
examples are parametric up-conversion, the generation of single microwave photons 
on demand [10, 12, 43], and squeezing of quantum states [44]. 

4. Conclusions 

We have fabricated superconducting three-junction flux qubits using electron beam 
lithography and two- angle shadow evaporation. The qubits are based on Al/AlO^/Al 
Josephson junctions with areas down to 0.02 fim 2 and high current densities above 
lkA/cm 2 . They show a level splitting above 4 GHz at 50 mK. We also fabricated 
readout dc SQUIDs and additional on-chip circuit elements such as shunting capacitors 
for the readout SQUID, transmission lines for the application of microwave signals to 
the qubits, and coplanar waveguide resonators. The latter are based on Nb films and 
have resonance frequencies ranging between about 1 and 10 GHz. They show internal 
quality factors above 10 5 at 1.5 K. 

We also studied the tunability of superconducting charge and flux qubits by 
external control fields. We derived the effective dipole moments of the qubits and 
show that they are by several orders of magnitude larger than those of natural atoms 
because of the much larger geometrical size of the artificial solid-state atoms (qubits). 
Furthermore, we derived the first- and second-order energy shifts of the qubit levels 
due to control fields as well as expressions for the matrix elements describing the 
level transitions at and away from the degeneracy point. We finally addressed some 
fundamental symmetry aspects of superconducting qubits. We argue that at the 
degeneracy point the qubit potential has mirror symmetry resulting in selection rules 
for level transitions equivalent to those in natural atoms. However, the symmetry of 
the qubit potential can be easily broken by external control fields. Then, the strict 
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selection rules are no longer valid. In particular, two-photon transitions are allowed 
resulting in an interesting up-conversion dynamics in a coupled qubit-resonator system, 
which we could observe in a recent experiment. 



Acknowledgments 

This work is supported by the German Science Foundation via SFB 631 and the 
German Excellence Initiative via the Nanosystems Initiative Munich (NIM). 

References 

[1] Nielsen MA and Chuang IL 2000 Quantum Computation and Quantum Information (Cambridge 
Univ. Press) 

[2] Deppe F, Mariantoni M, Menzel EP, Marx A, Saito S, Kakuyanagi K, Tanaka H, Meno T, Semba 
K, Takayanagi H, Solano E and Gross R 2008 Nature Physics 4 686 

[3] Kakuyanagi K, Meno T, Saito S, Nakano H, Semba K, Takayanagi H, Deppe F and Shnirman 
A 2000 Phys. Rev. Lett. 98 047004 

[4] Deppe F, Mariantoni M, Menzel EP, Saito S, Kakuyanagi K, Tanaka H, Meno T, Semba K, 
Takayanagi H and Gross R 2007 Phys. Rev. B 76 214503 

[5] Clarke J and Wilhelm FK 2008 Nature 453 1031 

[6] Schoelkopf RJ and Girvin SM 2008 Nature 451 664 

[7] Wallraff A. et al. 2004 Nature 431 162 

[8] Chiorescu I. et al. 2004 Nature 431, 159 

[9] Blais A, Huang RS, Wallraff A, Girvin SM and Schoelkopf RJ 2004 Phys. Rev. A 69 062320 
[10] Mariantoni M et al. 2005 eprint arXiv:cond-mat/ 0509737 (unpublished). 

[11] Johansson J, Saito S, Meno T, Nakano H, Ueda M, Semba K and Takayanagi H 2006 Phys. Rev. 

Lett. 96 127006 
[12] Houck AA et al. 2007 Nature 449 328 

[13] Sillanpaa MA, Park JI and Simmonds RW 2007 Nature 449 438 
[14] Majer J et al. 2007 Nature 449 443 
[15] Astanev O et al. 2008 Nature 449 588 
[16] Hofheinz M et al. 2008 Nature 454 310 
[17] Fink JM et al. 2008 Nature 454 315 

[18] Mariantoni M, Deppe F, Marx A, Gross R, Wilhelm FK and Solano E 2008 Phys. Rev. B 78 
104508 

[19] Mabuchi H and Doherty AC 2002 Science 298 1372 

[20] Haroche S and Raimond JM 2006 Exploring the Quantum: Atoms, Cavities, and Photons 
(Oxford Univ. Press) 

[21] Walther H, Varcoe BTH, Englert BG and Becker T 2006 Rep. Prog. Phys. 69 1325 
[22] Cohen- Tannoudji C, Diu B and Laloe F 1977 Quantum Mechanics (Wileylnterscience, New 
York) 

[23] Mooij JE, Orlando TP, Levitov LS, Tian L, van der Wal CH and Lloyd S 1999 Science 285, 
1036 

[24] Orlando TP, Mooij JE, Tian L, van der Wal CH, Levitov LS, Lloyd S and Mazo JJ 1999 Phys. 
Rev. B 60 15398 

[25] van der Wal CH, ter Haar ACJ, Wilhelm FK, Schouten RN, Harmans CJPM, Orlando TP, Lloyd 

S and Mooij JE 2000 Science 290 773 
[26] Dolan GJ 1977 Appl. Phys. Lett. 31 337 

[27] Ambegaokar V and Baratoff A 1963 Phys. Rev. Lett. 10 486 

[28] Deppe F, Saito S, Tanaka H and Takayanagi H 2004 J. Appl. Phys. 95 2607 

[29] Wallraff A et al. 2007 Phys. Rev. Lett. 99 050501 

[30] Frunzio L, Wallraff A, Schuster D, Majer J, Schoelkopf R 2005 IEEE Trans. Appl. Supercond. 
15, 860 

[31] Pozar DM Microwave Engineering 2005 (John Wiley Sz Sons) 
[32] Girfoni M, Hanggi P 1998 Phys. Rep. 304 229 

[33] Meschede D, Walther H. and Muller G 1985 Phys. Rev. Lett. 54 551 
[34] Rempe G, Walther H. and Klein N 1985 Phys. Rev. Lett. 58 353. 
[35] Raimond JM, Brune M. and Haroche S 2001 Rev. Mod Phys. 73 565 



Superconducting Quantum Circuits 



20 



[36] Yoshihara F, Harrabi K, Niskanen AO, Nakamura Y and Tsai JS 2006 Phys. Rev. Lett. 97 
167001 

[37] Averin DV and Bruder C 2003 Phys. Rev. Lett. 91 057003 

[38] Sillanpaa MA, Lehtinen T, Paila A, Makhlin Yu, Roschier L and Hakonen PJ 2005 Phys. Rev. 
Lett. 95 206806 

[39] Duty T, Johansson G, Bladh K, Gunnarsson D, Wilson C and Delsing P 2005 Phys. Rev. Lett. 
95 206807 

[40] Konemann J, Zangerle H, Mackrodt B, Dolata R and Zorin AB 2007, Phys. Rev. B 76 134507 
[41] Liu Y-X, You JQ, Wei LF, Sun CP and Nori F 2005 Phys. Rev. Lett. 95 087001 
[42] Brune M, Raimond JM, Haroche S 1987 Phys. Rev. A 35 154 
[43] Liu Y-X, Wei LF and Nori F. 2004 Europhys. Lett. 67, 941 
[44] Moon K and Girvin SM 2005 Phys. Rev. Lett. 95 140504 



